The Effect of Massive Perturbers on Extreme Mass-Ratio Inspiral Waveforms 
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Extreme mass ratio inspirals, in which a stellar-mass object orbits a supermassive black hole, 
are prime sources for space-based gravitational wave detectors because they will facilitate tests of 
strong gravity and probe the spacetime around rotating compact objects. In the last few years of 
such inspirals, the total phase is in the millions of radians and details of the waveforms are sensitive 
to small perturbations. We show that one potentially detectable perturbation is the presence of 
a second supermassive black hole within a few tenths of a parsec. The acceleration produced by 
the perturber on the extreme mass-ratio system produces a steady drift that causes the waveform 
to deviate systematically from that of an isolated system. If the perturber is a few tenths of a 
parsec from the extreme-mass ratio system (plausible in as many as a few percent of cases) higher 
derivatives of motion might also be detectable. In that case, the mass and distance of the perturber 
can be derived independently, which would allow a new probe of merger dynamics. 



I. INTRODUCTION 

Space-based gravitational wave detectors such as the 
Laser Interferometer Space Antenna (LISA) are expected 
to see a wide variety of sources in their ~ 10 -4 — 1CP 1 Hz 
sensitivity band. Of these, extreme mass-ratio inspirals 
(EMRIs), a stellar-mass compact object (SCO) spiraling 
into a supermassive black hole (SMBH), are considered 
particularly promising because they can probe strong 
gravity over millions of radians of phase evolution in the 
last few years of evolution. As a result, EMRI waveforms 
serve as highly precise probes of strong gravity and of the 
spacetime around rotating SMBHs. Considerable study 
has been devoted to astrophysical scenarios for EMRIs 
[IHH as well as to the analysis of their waveforms (34l0| . 

There has been less exploration of the possibility of 
deviations from isolated EMRI waveforms that might oc- 
cur due to environmental effects (see e.g. 0, HTI - fl4| for a 
study of differences caused by an accretion disk around 
the SMBH) . Here we point out an effect that has not been 
considered in this context: the acceleration of the EMRI 
system by a nearby (distance of roughly a few tenths of a 
parsec or less) secondary SMBH. As we demonstrate, this 
acceleration leads to phase drifts of fractions of a radian 
over a year of inspiral, which is potentially detectable 
from EMRIs of plausible signal strength. Depending on 
the fraction of galaxies that merge, and on the fraction of 
time in such mergers that the secondary SMBH is within 
a few tenths of a parsec of the primary, this could affect 
as many as a few percent of EMRIs. 

The detection of such an effect could yield a new probe 
of galactic merger dynamics, providing a measure of the 
ratio of the secondary SMBH's mass and its distance to 
the EMRI. If such effects arc not present in a detected 
gravitational wave (GW), then one can place an upper 



limit on the density of SMBHs inside some radius of a 
few tenths of a parsec. If this is the case, then one would 
confirm that, as far as LISA is concerned, EMRIs occur 
in vacuum. 

This paper is organized as follows: In Sec. |H] we do 
a simple analysis of the acceleration effect as it would 
apply to a signal of constant frequency and amplitude, 
which wc expand on in the Appendix. In Sec. IIIII we ex- 
plain how to model real EMRI waveforms, for the partic- 
ular case of quasi-circular, equatorial orbits, and explain 
how to implement modifications to model an accelera- 
tion effect. In Sec. IIVI we extend the simple analysis of 
Sec. |TT] to real waveforms and perform a dephasing and 
an overlap study. In Sec. |V] we explore whether some 
of these deviations can be masked by adjustments of 
EMRI system parameters. We present our conclusions 
in Sec. IVII In most of this paper, we use geometric units 
with G = c = 1. For reference, in this system of units, 
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II. SIMPLE MODEL 

Here we present the basic effects of acceleration in a 
simplified model. Wc assume that there is an EMRI of a 
SCO into a (primary) SMBH with mass M, on the x-y 
plane, with orbital and spin angular momentum in the z 
direction. We further simplify the scenario by assuming 
GWs of constant frequency and amplitude. Let us also 
assume there is a secondary SMBH in a circular orbit 
about the EMRI's center of mass (COM). Suppose that 
the secondary SMBH has a mass M Sec and the total mass 
of the system M Tot = M, + M Soc . Suppose also that 
the semi-major axis of the circular orbit of the primary- 
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FIG. 1. Schematic view of the EMRI system (in the xy plane) , 
the massive perturber M Soc (at a distance r So c from the EMRI 
SMBH), and the line of sight, i is the inclination between 
the primary-secondary SMBH's orbital angular momentum 
vector and the line of sight. 



secondary SMBH system is r Soc , and that it is inclined to 
the line of sight at an angle i (here i is zero for a face-on 
binary and 90° for an edge-on binary). We depict this 
scenario in Fig. [1] 

If these systems are well-separated, then the EMRI's 
COM will move essentially at a constant velocity relative 
to us, with a projection into our line of sight of 
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where i> Ncwt = {GM Tot / r^) 1 / 2 is the Newtonian virial 
velocity, cj Ncwt = (GM^t/r^J 1 ' 2 is the Newtonian an- 
gular velocity for an object in a circular orbit, and 5 is 
an initial phase, with (w Ncwt i + S) the orbital phase of 
the EMRI-Secondary system. A constant relative speed 
is entirely absorbed in a redefinition of the masses. As 
such, constant relative velocities cannot enter any of our 
results. 

If the EMRI's COM is sufficiently close to the sec- 
ondary SMBH, then the former will experience an mea- 
surable acceleration, which will produce a net Dopplcr 
phase drift relative to the best-fit waveform. Note that 
the orbital period for M Tot = 10 6-7 Mq and r Soc = 
0.1 — 1 pc is at least ~ 10 3-5 years, so for the duration of 
a LISA observation the binary will not change phase sig- 
nificantly. Tidal effects on the EMRI system due to the 
perturber can be neglected, as this acceleration scales as 
the inverse cube of r So( . (see Sec. lIIIBl for more details). 

If the EMRI system is accelerated by an amount v 
relative to its original line of sight speed over a time t, 
then the GW phase difference compared to the initial 
frequency is A<I> GW = ^vtN/c, where N is the number 
of radians in the waveform (see Appendix [A] for a more 
detailed explanation of this effect). Let us designate by e 
the detectable fractional phase shift: e = A3> GWj detect /N. 



As a fiducial value we will use e = 10~ 7 , or 0.1 radians 
over ~ 10 6 radians for a typical one-year inspiral. To 
leading-order in a Taylor-expansion about w Ncwt t = 0, we 
then have 

1 i>t I , . w . r . M Soc t 

-— = - sini (sm<5) =e, (2) 

2 c 2 r Scc r Soc 

where we have here neglected a constant term that is 
non-observable. Solving for the distance at which this is 
satisfied gives 
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The next order term in the phase shift scales as 
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Setting this equal to e and solving for r, we find 
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Additional corrections can be computed similarly. 

Therefore, for BH masses <; 10 6 Mq and separations 
of a few tenths of a parsec or less, acceleration can cause 
a detectable shift in the simplified waveform. As we find 
in Sec. IIV| this shift is proportional to the combination 
A = M Sec /r| isc . For separations of a few hundredths of 
a parsec or less, higher order derivatives are measurable. 
In this case, the detectable shift in the waveform is cap- 
tured by the linear combination of A and other highcr- 
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order derivative terms, such as B 
a sufficiently small r Scc one could then measure both A 
and B and thus disentangle M Sec from r Soc . 

The range of masses and separations that could be ob- 
served, given a sufficiently strong EMRI-pcrturber sys- 
tem are depicted in Fig. [5] In this figure, we show with 
solid lines the constraint given by Eq. @, and with a 
dashed line that of Eq. (with M. = M Seo for sim- 
plicity), where the black, red and blue colors correspond 
to e = 10~ 7 , 10~ 6 and 10~ 5 . A larger value of e corre- 
sponds to more conservative choices of what is detectable 
by LISA. The area above the curves show the values of 
M Soc and r Soc that could be measurable. For compari- 
son, we also show the region of (M Soc ,r Soc ) space that 
fall in the pulsar-timing-array (PTA) sensitivity band. 
Of course, for PTAs to individually resolve such bina- 
ries, their distance to Earth would have to be sufficiently 
small 15 1 . In principle, however, this scenario allows for 
the possibility of coincident future detection of GWs with 
LISA and PTAs. 
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FIG. 2. Range of secondary masses and separations that could 
be measurable by LISA given a sufficiently strong EMRI. The 
region above the solid and dashed lines would be observable. 
Measurement of the leading-order effect gives a determina- 
tion of the combination (M Sec sin t) /r| BC , while measuring the 
next-order effect gives a determination of the combination 



(M Soc sini) ' / 



3/2 / 7/2 



Thus measuring both effects together al- 



lows both M Soc sin i and r Sec to be determined. 



We now discuss the detectability of the changes dis- 
cussed above for a realistic EMRI. In Sec. I VII we return 
to the question of how common it will be to have a sec- 
ondary SMBH this close. 



III. REALISTIC EMRI WAVEFORMS 
A. Standard EOB Modeling 

We employ the effective-one-body (EOB) formalism to 
model waveforms with and without the acceleration cor- 
rection. This formalism was initially developed in [HI, [13] 
to model comparable-mass BH binary coalescences. Im- 
provements and extensions to other binaries were de- 
veloped in [l8r]29| and compared to a set of numerical 
relativity results in [3C - 32 1 and to self- force calculations 
in [33|, Hi| . Recently, |E - 1C | combined the EOB approach 
with BH perturbation theory results to model EMRI 
waveforms for LISA data-analysis. We here concentrate 
on the formulation of [l(|, as it is applicable to EMRIs, 
the systems of interest in this paper. 

We focus on quasi-circular EMRI inspirals in the equa- 
torial plane of a spinning BH because they are simpler to 
model. We define the following orbital parameters: the 
SCO's mass to*; the SMBH's mass M.; the total mass 
M = m* + M.; the reduced mass \i = m+m,/M] and the 
symmetric mass-ratio v = u/M. We also assume that the 
EMRFs orbital angular momentum is aligned with the 
MBH's spin angular momentum S, = a. M, = q,M?, 
where a. = S m /M m is the MBH's spin parameter and 
<7» = o, % jM % is its dimcnsionlcss spin parameter. We em- 
ploy the adiabatic approximation, in which we assume 



that the radiation-reaction time-scale is much longer than 
the orbital one. 

With this at hand, let us now describe the EOB ap- 
proach we employ. In the adiabatic approximation, the 
GW phase can be obtained by solving 



<j) = UJ, 



(6) 
(7) 



where ui = </> is the orbital angular frequency, with <j) the 
orbital phase, overhead dots stand for time derivatives, 
E is the system's total energy and F is the GW energy 
flux. The energy of the system is [HI 



E = M. + m* 



1 - 2M./r±q. AI? /2 /r 3 / 2 
1 - 3M./r ± 2q. M? /2 /r 3 / 2 



(8) 



where the ± stands for progradc or retrograde orbits. In 
this equation and all throughout the rest of this paper, 
we ignore sub-leading corrections that are proportional 
to the EMRFs mass-ratio. In practice, this means we ig- 
nore conservative and second-order dissipativc self-force 
effects, i.e. the effect of the SCO on its own geometry, as 
well as the SCO's spin. 

The GW flux can be written in the factorized form 
of OH H3 HI HI , which in the adiabatic regime is 



1 8 1 
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where the multipolc-decomposed waveforms are 

w«) = h Z t,ep S%n T ^ e lS ^ ( Pim ) e , 



(9) 



(10) 



and where e p is the parity of the waveform (i.e., e v = 
if t + to is even, e p — 1 if £ + in is odd). The quantities 
(5^(w), Ttwjv), 5irn.(v) and pe m {v)) in Eq. (jTOj) can be 
found in 24, 53, r2al • The Newtonian waveform is 




nfe ) c l+tp v e+t * Y t _ 



.(tt/2,0). (11) 



where 



are spherical harmonic functions, while 
( and Q +ep are numerical coefficients [2? 



We enhance the flux of Eq. © by linearly adding BH 
absorption terms and calibration coefficients that are fit- 
ted to a more accurate, numerical flux jTjjj]. The first 
modification is necessary as BHs lose energy due to GWs 
that both escape to infinity and fall into BHs. The second 
modification accounts for the fact that the bare fluxes 
written above are built from low-velocity (PN) expan- 
sions, and as such, are not sufficiently accurate by them- 
selves for long evolutions, even after the resummations 
introduced. 

The above differential system is solved with the post- 
circular initial conditions of II 71 . enhanced with a mock- 
evolution at 100M, (see e. g. PI|). The orbital phase can 
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then be used in the waveforms of Eq. (fTUf , together with 
the fact that for quasi-circular orbits 



[1 - gl (M.cu)} 2/3 
(M.w) 2 / 3 



(12) 



where r is the EMRI's separation and v = (M.w) 1 / 3 by 
Kepler's third law. With the waves at hand, we then 
compute the GW phase and its amplitude via 
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(13) 



The GW phase as defined above needs to be unwrapped 
every 27T, so in practice it is simpler to define the time 
derivative of this quantity and then obtain via inte- 
gration. 



scales away from the COM, and one parsec translates to 
~ 10 13 M Q in geometric units. Since the tidal force scales 
as the inverse of the separation cubed, any tidal effects 
are insignificant. 

Given that this type of generation effects are sup- 
pressed, are there any others that should be included? 
The dominant generation effect is simply a Dopplcr shift 
in the frequencies, which then leads to an integrated mod- 
ification in the GW phase (see the Appendix for a de- 
tailed explanation of this Doppler effect). In this sense, 
such a correction is similar to the integrated Sachs- Wolfe 
effect for GWs [37|, where here the perturbation to the 
potential is given by a third body, instead of some cos- 
mological background. The implementation of this cor- 
rection to an EOB evolution is simple: divide the right- 
hand-side of Eq. (|7|) by the appropriate Doppler factor 



B. Modifications to EOB Modeling 

How do we incorporate the effects of an external accel- 
eration into GW modeling within the EOB framework? 
Let us first distinguish between wave generation and wave 
propagation effects. By the former, we mean effects that 
arise in the near-zone (less than a gravitational wave- 
length away from the EMRI's COM) and that generate 
GWs due to the inspiral of the EMRI. By the latter, we 
mean effects that arise after the system has generated a 
GW and it propagates out to the wave- zone, where the 
observer is located, many gravitational wavelengths away 
from the source. 

As is expected, all propagation effects, such as the 
backscattcring (or tails) of GWs off the metric of the 
secondary SMBH, occur beyond Newtonian (leading) or- 
der in post-Newtonian theory [36| , and can be safely ne- 
glected here. The presence of an external source, how- 
ever, does introduce non-negligible modifications to the 
generation of GWs. One could incorporate such effects 
by introducing an external, vectorial force to Hamilton's 
equations in the direction of the perturber. This force 
would simply be the product of the the total mass of the 
system and the time derivative of the velocity of Eq. (fTJ . 
The modeling of this effect would require a non-adiabatic 
evolution, i.e. the evolution of the full set of Hamilton's 
equations, without assuming circular orbits or using Ke- 
pler's third law. One expects that such force would in- 
duce eccentricity and inclination in the inspiral, driving 
the SCO out of the equatorial plane of the secondary 
SMBH. 

One can estimate the magnitude of this effect by con- 
sidering the tidal force effect of the perturber on the 
COM relative to the SCO's acceleration due to the sec- 
ondary SMBH. Since the tidal force scales as F Tidal = 
M Soc /V 3 oc , this effect is suppressed relative to the ac- 
celeration by a factor of r/r Scc ~ 10~ 4 for an EMRI 
with orbital separation of 30to, and a primary-secondary 
SMBH orbital separation of 0.01 pc. The ratio is this 
small because the perturber is assumed to be at parsec 



4> = uj -> = oj [1 + v los (t,S = tt/2)] . (14) 

In this equation, we have not included the appropriate 
Lorentz factor T, since t> Nowt /c -C 1, and we can linearize 
in this quantity. Moreover, we have removed the constant 
velocity drift component of u los by choosing S = n/2, as 
the former is not measurable. 



IV. PERTURBING ACCELERATION EFFECT 
ON RELATIVISTIC EMRI WAVEFORMS 

A. Preliminary Considerations 

With the machinery described in Sec. IIII1 we can con- 
struct modified EMRI waveforms as a function of time, 
for a given value of the second MBH mass and separa- 
tion to the EMRI's COM. We consider the following two 
EMRI systems, integrated for one year each: 

• System I: The primary SMBH has mass m, = 

10 5 Mq and spin parameter q, = 0.9, while the 
SCO has mass and spin parameter = 10 Mq 
and q± = 0. This system inspirals for ^ 6 x 10 5 
rads of orbital phase between orbital separations 
r/M £ (16,26). In this range the orbital veloci- 
ties are v £ (0.2, 0.25) and the GW frequencies are 
/ow G (0.005, 0.01) Hz. 

• System II: The primary SMBH has mass m, = 

10 6 Mq and spin parameter <j. = 0.9, while the 
SCO has mass and spin parameter m* = 10 Mq 
and = 0. This system inspirals for ~ 3 x 10 5 rads 
of orbital phase between orbital separations r/M £ 
(11, r isco)- In this range the orbital velocities are 
v £ (0.3, fisco) and the GW frequencies are / Gw £ 
(0.001,/™°) Hz. 

System I exits the most sensitive part of the LISA band at 
around 16M, which is why we stop the evolution there. 
In contrast, Sys. II is stopped when the SCO reaches 
the innermost stable circular orbit (ISCO). For each of 
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TABLE I. Summary of System properties. All masses are in 
units of M© , r So c is in units of parsecs and all time scales are 
in units of years. The time to merger is here estimated as 
T G w = r/ f GW , where r GW is the rate of change of the orbital 
separation due to GW emission and for r we take the values 
in the itemized list. The superscript star stands for quantities 
associated with the SCO-SMBH system, while the solid dot 
stands for those associated with the SMBH-SMBH system. 



these systems, we explore a variety of secondary SMBH 
masses M Soc = (10 5 , 10 6 )M Q as well as a variety of sepa- 
rations r Soc = (0.01,0.1,1) pc. Larger secondary masses 
are also possible; these would have equivalent effects on 
the EMRI at correspondingly larger distances r ~ M 1 / 2 . 
(For example, M Sec = 10 9 M Q at r — 30 pc would have 
equivalent effects to M Seo = 10 6 M Q at r = 1 pc.) All of 
this information is summarized in Table HI including the 
orbital periods T 0rh and the time to merger due to GW 
emission T GW . In all cases we set sint = 1 and S = n/2, 
as this leads to the largest possible effect. The reasoning 
behind this is that if this effect is not observable with 
this choice of parameters, it will not be observable with 
any other choice. 
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FIG. 3. Dephasing for Sys. I and II as a function of time in 
units of months, for a variety of separations and masses of the 
perturber. 



pc] , while different line styles correspond to different per- 
turber masses [M Soc = (10 5 , 10 6 )M Q ]. Observe that for 
both Systems a dephasing of order 0.1 rads is achieved 
for separations r Scc < 0.1 pc over less than one year. 
This is consistent with the estimates of Sec. HU Sim- 
ilarly, more massive perturbers enhance the dephasing 
roughly by one order of magnitude. The amplitudes of 
the waveforms are not shown in this figure; they disagree 
at the level of 10~ 3 for Sys. I and 10~ 4 for Sys. II. 

The magnitude and shape of the dephasing depends 
on how far away and massive the perturber is. One can 
show that the dephasing scales as A$ GW oc N M Soc T /r 2 cc , 
where T is the observation time and is the number of 
cycles. Since there is a factor of two fewer GW cycles 
in Sys. II relative to Sys. I, then the dephasing for the 
former is also smaller by a factor of two. 



B. Dephasing Study 

Let us define the dephasing between waveforms as fol- 
lows: 



(f>»°- Acc 



(15) 



where is the GW phase of an EMRI waveform with 
an accelerated COM, while $ G ° w Acc is that of an iner- 
tial COM. We have here aligned the waveforms in time 
and phase before computing this dephasing. This align- 
ment is equivalent to minimizing the statistic in Eq. (23) 
of [3l| , which in turn is the same as maximizing the fitting 
factor over time and phase of coalescence in a matched 
filtering calculation with white noise 3 ill . The alignment 
is done here in the same way as in [8 4l0| . 

Figure [3] plots the dephasing of the dominant (£, m) = 
(2, 2) GW mode as a function of time in months for Sys. I 
and II. The different line colors/shades correspond to dif- 
ferent separations to the perturber [r Soc = (0.01,0.1,1) 



C. Overlap Study 

The dephasing study of the previous subsection is 
suggestive, but not sufficiently quantitative to assess 
whether such types of corrections are large enough to 
be measurable. Let us then perform a slightly more so- 
phisticated data analysis study here. 

Given any time series A(t) and B(t), one can construct 
the inner-product 



(A\ B) = 4Re / 
Jo 



MS) B*(f) 

Sn(f) 



df 



(16) 



where the overhead tildes stand for the Fourier trans- 
form, the star stands for complex conjugation and S n (f) 
is the spectral density of noise in the detector. We choose 
here the sky-averaged version of the noise curve presented 

in mm. 

With this inner-product, we can now construct some 
data analysis measures. The signal to noise ratio (SNR) 
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of signal A is 



p(A) = ^/{AjAj , 
while the overlap between signals a and b is 

(A\ B) 



M = max- 



A) (B\ B) 



(17) 



(18) 



with the mismatch MM = 1 — M. The max label in 
Eq. (|18[) is to remind us that this statistic must be maxi- 
mized over an event time (e.g., the time of coalescence of 
the EMRI system) and a phase shift 126] . If the overlap 
is larger than 97% (or equivalently, if the mismatch is 
lower than 3%), then the difference between waveforms 
A and B is sufficiently small to not matter for detection 
purposes (see e. g. |40|). The minimum overlap quoted 
above (97%) is mostly conventional, motivated by the 
fact that the event rate scales as the cube of the overlap 
for a reasonable source distribution. For an overlap larger 
than 97%, no more than 10% of events are expected to 
be lost at SNRs of O(10). Of course, for larger SNRs, 
one might not need such high overlaps, although EMRI 
sources are expected to have SNRs < 100. 

Whether the difference between waveforms A and B 
can be detected in parameter estimation can be assessed 
by computing the SNR of the difference in the waveforms 
Sh = A - B: 



p{5h) = sf(Jh\ Sh) = 4 Re 



Sh(f) Sh (f) 

Sn(f) 



df. (19) 



When this SNR equals unity, then one can claim that A 
and B arc sufficiently dissimilar that they can be differ- 
entiated via matched filtering (see e. g. |41|). 

We applied these measures to EOB waveforms with 
and without acceleration of the COM. The results arc 
plotted in Fig. H] function of observation time in 
months. The vertical dotted lines correspond to obser- 
vation times of (0.5,2,4,6,9,12) months, and the num- 
bers next to them, in parenthesis, stand for the SNR 
of Sys. I and II for that observation time. The differ- 
ent line styles and colors correspond to mismatches and 
SNRs of the error for different secondary systems. Ob- 
serve that the mismatch is always smaller than 0.03 (the 
solid black horizontal line), suggesting that this effect 
will not affect detection. Observe also that the SNR of 
the difference reaches unity (the dashed black horizontal 
line) in between 6 and 12 months of observation, and for 
the Msec = 10 6 M Q , p(<5h) reaches ~ 10 after one year. 
This suggests that given a sufficiently strong EMRI with 
SNR ~ 50 — 100, the magnitude of this effect is in prin- 
ciple detectable within one year of coherent integration. 



V. DEGENERACIES 

Now that we have determined that there exists a set 
of plausible pcrturbcr parameters for which the magni- 
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FIG. 4. Mismatch as a function of time in units of months 
for Sys. I and II and different perturber masses, all at a sep- 
aration of r So c =0.1 pc. SNRs for Sys. I and II are given in 
parentheses for a source at 1 Gpc. 



tude of the correction could be measurable, let us con- 
sider the possibility of degeneracies. That is, let us in- 
vestigate whether we can mimic an acceleration of the 
COM by changing the intrinsic parameters (the com- 
ponent masses, the spin parameter, etc.) in the non- 
accelerating waveform. The simplest way to see whether 
this is possible is to study the frequency dependence of 
the GW modification introduced by the COM's acceler- 
ation. 

Let us then remind ourselves of how the frequency- 
domain representation is constructed. For this, we em- 
ploy the stationary- phase approximation (see e. g. [42j). 
under which, the frequency-domain waveform is simply 



Hf) = Af 



-7/6 e ^(f) 



(20) 



where the Newtonian (leading-order) amplitude is A — 
^-2/330-1/2 m 5/6 D1 1 , with M = ?y 3 / 5 M, while the 
phase is constructed from 

lK/)=~+27r/t(/)-20(/), (21) 

where the second term arises due to the Fourier transform 
and the third term due to the oscillatory nature of the 
time-domain waveform. 

The phase of the frequency-domain waveform in the 
stationary phase approximation is then controlled by 
these last two terms in Eq. (|2"Tj) . The first term can be 
computed via 



2tt/*(/) = 2tt/ 



//2 



r(F') 
F' 



dF' . 



(22) 



where / is the GW frequency, while the second term can 
be calculated from 



<Kf) = ^ 



in 



r(F')dF' , 



(23) 
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where t(F) = F/F and F is the orbital frequency. 

The Doppler correction to the waveform comes in the 
calculation of <fi(f), as this is simply the integral of the 
frequency. For simplicity, we can reparameterize the 
Vioe(t) v 1ob (F), by noting that, to Newtonian order, 



27rF(t) 



4-3/2 
M 



— (t c - t) 

5M v ; 



which we can invert to obtain 

, . 5M . 



-3/8 



•8/3 



(24) 



(25) 



where t c is the time of coalescence of the EMRI system. 
Taylor-expanding Eq. ([1} about w Nowt i = 0, we find 



v lOB ~v + v! (2nMF) 8/3 + v 2 (2nMF) 16/3 



v 3 {2itMF) 



-s 



(26) 



where the vi coefficient are the following frequency- 
independent functions: 



I'D 



1/2 / M Seo x 1//2 



1 M Tot t 2 c 1 /M T 

2^ C ° S<5+ 2 



cos 5 — 



Mrj, 



1/2 



sin<5 



3/2 t . 



sin 8 



10 MM S , 
512 r 2 



-77 (sin i) 



r 2 



lM Tot t 2 c . ' 
-- — = — -sins 

2 r| ec 

25 M 2 M S(!C 
131072 HL 



sin <5 



Mr 



tec 



sin l 



\ r Se c 



1/2 



■ cos S 



1/2 



cos -I — sm 



V3 



xr] sm 1 , 

125 M 3 M 8ec M T 



100663296 



-i] sinismt. 



(27) 



Notice that v is of 0(M 1 / 2 /r^J), v x is of 0(M 2 /rjJ 7 

v 2 is of 0{M 7 ' 2 /rl'J) and v 3 is of C(Af 5 /^Lc)- 

With these relations at hand, we can now compute the 
correction to the frequency-domain waveform phase in 
the stationary phase approximation. Denoting by Atp = 

V'Acc-V'no. Ace we find that 

rf/2 

A^ = -4ir t(F') v }oe {F') dF' , 
5ttM r s/2 



24 77 



(2iTMF'y 8/3 



v + v 1 (2ttMF') 8/3 



v 2 (2ttMF') 16 + v 3 (2irMF') 8 



dF' 



(28) 



where in the second line we have used that to Newtonian 
order 



F 



48 



5nM' 



{2irMF) 



11/3 



(29) 



Normalizing this phase correction by the Newtonian form 
of the frequency-domain waveform phase, we find 



128 



(7rMf)- 5/3 



Setting S = tt/2 = t, the above expression simplifies to 



A^= JLfrMf)- 6/a 



25 MM, 



M SBC t c 4M SBC M Tot t 2 c 



3 rl c 9 ? 
25 MM Scc M Tot t 2 c 



1248 r 2 cc 2496 
125 M 2 M Soc M Tot t c 



(nMf) 



-8/3 



1032192 
625 



M 3 M,M T 



(nMf) 



-16/3 



1094713344 



(irMfY 



(31) 



Let us now discuss this result in more detail. The first 
two terms inside the square bracket in Eq. pip arise due 
to a constant misalignment between the time of coales- 
cence of the EMRI system and the primary-secondary 
SMBH system (we have implicitly set the latter to zero). 
This effect can be absorbed via a redefinition of the chirp 
mass, and thus, it is not observable. All other lines 
in Eq. pip , on the other hand, contain a non-trivial 
frequency-dependence and they cannot be reabsorbed via 
a redefinition of intrinsic parameters. 

A physical way to think about this is the following. 
Given a signal and a template without modeling a sec- 
ondary perturber, one would like to maximize the phase 
coherence by shifting the template's phase and frequency. 
Such a shift corresponds to an adjustment of the total 
mass and the chirp mass, which eliminates the first term 
in Eq. p"Tj) . Once this shift is done, however, there are 
no other template parameters that can be shifted, while 
the frequency derivatives of the signal and template will 
continue to disagree. 

The many terms that arise in the second, third and 
fourth lines of Eq. pjj are due to the Taylor expansion 
in w Ncwt rj, which in frequency space has become an expan- 
sion in inverse powers of (M /) and (r Scc /M Soc ). Clearly, 
the second line Eq. p"Tj) is dominant over all others as 
it scales with r~ 2 to leading order, while the third and 
fourth lines scale as r~^. If M Scc /r Soc is large enough, 
however, one might be able to measure the coefficients in 
front of both the dominant J -8 / 3 term and the J" 16 / 3 or 
/~ 8 term. This would then imply that one could break 
the degeneracy between M Soc and r Soc in the leading order 
term and measure both quantities. 

Ignoring such possible acceleration effects could intro- 
duce a bias in the extraction of parameters via matched 
filtering [4^, |44 |. Imagine, for example, that an EMRI 
GW is detected with vacuum templates. One would then 
proceed to extract parameters from this detection, such 
as the primary SMBH's and SCO's mass with some error 
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bars. Usually, the error estimate accounts for statisti- 
cal error plus possible systematics with the modeling. 
The acceleration effect here would be one such system- 
atic, whose magnitude would have to be determined via 
a careful Markov- Chain Monte-Carlo exploration of the 
likelihood surface. 

Notice also that in Eq. ([31]) we have kept only the New- 
tonian contribution to an infinite post-Newtonian expan- 
sion. This is essentially because in Eqs. (f2"4")l and (|29|) wc 
have dropped all but the leading order, Newtonian term. 
Interestingly, the correction terms that arise at leading 
order are dominant over the Newtonian piece, as they 
depend on high inverse powers of frequency (in partic- 
ular, higher than 5/3). This implies that if a detailed 
parameter estimation study were to be carried out, these 
post-Newtonian terms should be taken into account, as 
they contribute at the same order as the Newtonian term 
in an inertial frame. 

The dependence of the correction in Eq. (|3"Tj) on dif- 
ferent powers of the frequency suggests that these terms 
are non-degencrate with the standard ones that appear 
in the non-accelerating GW phase. More precisely, the 
GW phase in an inertial frame is given by (see eg. [42|). 



2nft c -cj) c + — (nMfy 5/3 



(32) 



1 



'3715 55 



V 756 9 



where <j> c is the phase of coalescence and the ellipses stand 
for higher order terms in the post-Newtonian series. No- 
tice that there are no powers of J -8 / 3 , /~ 16 ^ 3 or /~ 8 in 
the above equation. Thus, the correction computed in 
Eq. (|3"Tj) is weakly correlated to the GW phase in an in- 
ertial frame, i. e. the off-diagonal elements of the Fisher 
matrix are small for the M Sec /r| oc coordinate sector rel- 
ative to the diagonal term. Although these results are 
suggestive, a more detailed analysis should be carried 
out to determine the level of correlation between all pa- 
rameters and the accuracy to which M Soc and r Scc could 
be extracted. 

Although the correction due to the acceleration of the 
COM seems to be weakly correlated to other intrinsic 
parameters, one might wonder whether it is degenerate 
with other effects not included in vacuum GR waveforms. 
Takahashi and Nakamura [45[ have studied the effect of 
the acceleration of the Universe in the frequency-domain 
form of the waveform. They find that 



Aip 



128 



(nMf)- 5/3 ^-Mz(irMf)- 8/S 



768 



(33) 



One can clearly see that this cosmological effect is de- 
generate with the one computed here [the second line 
in Eq. (|3ip ]. However, the magnitude of Eq. (f3"3"f is 
much smaller than that of Eq. (f3Tj). simply because 
H <C Af Scc /r| cc for all relevant perturbers considered 
here. For example, at small redshift, Hq ~ 10~ 23 km -1 
in geometric units, while at r Sec = 0.1 pc and for a 10 6 Af Q 



perturber, M Scc /r 2 oc ~ 10 19 km . The perturber sep- 
aration at which these effects become comparable is ap- 
proximately r Soc ~ 11 pc [Msec/^Mg)] 1 / 2 . 

Another possible source of degeneracy could be if 
there are corrections to general relativity that induce 
phase modifications with the specific frequency depen- 
dence found in Eq. ([31]) . In fact, we see that the result 
obtained here can be mapped to the parameterized post- 
Einsteinian framework \4M with the choice 



a = 0, 



25 M Scc M 
1248 r^r : 



3' 



(34) 



to leading order in M Soc /r Soc (see e.g. Eq. 1 in j44|). As 
found in that paper, however, there are no known alter- 
native theories to date that could potentially lead to the 
frequency dependence found in Eq. (|31l) . 



VI. DISCUSSION AND CONCLUSIONS 

We have shown that a ~ 10 6 M Q secondary SMBH 
within a few tenths of a parsec of the EMM system 
can produce detectable modifications in the waveform. 
A more massive secondary SMBH at a correspondingly 
larger distance would produce equivalent effects. It is not 
possible to say with certainty how common this will be. 
A rough upper limit can be obtained from the following 
observation. Since a redshift of z = 1 (corresponding 
roughly to 10 10 years), tens of percent of Milky Way-like 
galaxies have had a major merger (4(| HtJ . If the typical 
merger takes hundreds of millions of years, then at most 
a few percent will be involved in a merger at any stage. 
The fraction of time spent at separations ^ 1 pc remains 
uncertain; although there are well-understood dynamical 
processes that can reduce the secondary SMBH's separa- 
tion to ~ 1 pc and gravitational radiation will bring the 
binary to merger from ~ 10~ 3 pc, the transition between 
the regimes is uncertain (this is commonly called the "fi- 
nal parsec problem"; see, e.g., [48[ for a discussion). It 
is therefore possible that the system spends considerable 
time at roughly the detectable separations. 

We also note that when a secondary SMBH comes 
within a few tenths of a parsec of the primary, various 
dynamical effects temporarily increase the rate of close 
encounters of stellar- mass objects with both SMBHs [49j . 
As a result, it may be that a disproportionate number of 
EMRIs occur with a secondary SMBH nearby Indeed, 
recently [5(J estimated that more than 10% of all tidal 
disruption events could originate in massive black hole 
binaries, so if the EMRI fraction is similar it corresponds 
to our rough estimate. 

In these cases, measurement of an EMRI phase shift 
affords a new way to detect the presence of a binary 
SMBH. If the separation is close enough to measure an 
additional derivative of the motion, then the degeneracy 
between the secondary mass and its distance is broken. If 
the EMRI-SMBH system is sufficiently close, then pulsar 
timing measurements fl5| might also be able to detect 



9 



gravitational waves from the SMBH-SMBH binary. Al- 
ternatively, if no phase shift is detected, then this implies 
that there are no secondary SMBH in a radius of a few 
tenths of a parsec, thus implying an upper limit on the 
density of BHs close to the detected EMRIs. In princi- 
ple, therefore, EMRIs have another astrophysical link in 
addition to their utility in testing general relativity. 

The importance of the astrophysical environment in 
EMRI GW modeling is a double-edged sword. Although 
on the one hand, one could potentially extract some as- 
trophysical information, on the other, these effects could 
make it difficult to test general relativity [H)]. For such 
tests to be possible, one must have complete control of 
the waveforms within general relativity. If the astrophys- 
ical environment needs to be included, then the modeling 
might be dramatically more difficult. We note here, how- 
ever, that only a fraction of EMRIs would experience the 
astrophysical environment effect discussed here. If devi- 
ations from general relativity are present, on the other 
hand, these should be present for all EMRIs. Thus, in 
principle, a statistical analysis would allow us to disen- 
tangle deviations in our waveforms to discern whether 
they have an astrophysical or theoretical origin. 
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Appendix A: Acceleration Effect 

Here we explain in more detail how the Doppler correc- 
tion to the waveform comes about. Let us consider the 
effect of an acceleration on the COM position vector. For 
simplicity we consider the toy-model of a perfect circular 
orbit with angular velocity cj, whose position vector in 
the COM can be parameterized as 



x = b (cos tut , sin uit, 0) , 



(Al) 



where b is the binary's separation and we have erected 
a Cartesian coordinate system, with the binary in the 
x-y plane. If an external force is present that causes 
an acceleration, this in turn will cause a displacement 
x — > x 1 = x + Sx. Let us parameterize the magnitude 
of this displacement as \Sx\ = (l/2)v loa t 2 , which holds to 
Newtonian order for a uniformly accelerated body. One 



can then show that the shift in the magnitude of the 
COM velocity vector is simply 



= W\ + \v lo »t (x ■ Sx) + 0{vlf) , (A2) 



where v = x is the unperturbed velocity vector. Notice 
that there is a factor of 1/2 here, just as in the estimates 

of scc.ini 

Before proceeding, it is useful to concentrate on this 
velocity shift. Choosing <5 = 7r/2 = l, one can easily 
show that 



Mr 



"%cwt (^NcwtO 



(A3) 

upon Taylor expanding about w Nowt i <C 1. We can take 
the time-derivative of v los and then Taylor-expand again 
to find: 



Mr„ t 



^Ncwt (^NcwtO 



(A4) 

Obviously, this is the same as simply Taylor-expanding 
v los to leading order. 

One effect of the COM velocity drift is a Doppler shift 
to the waveform. Special relativity predicts that if a fre- 
quency source is moving with velocity v away from the 
observer at an angle 9, then the frequency observed is 



J = ^ (1 + vcos6) 1 



1 — v cos 9 + v [ cos 



(A5) 



where uj is the frequency of the source, u/ is the fre- 
quency the observer detects and T = (1 — u 2 ) -1 / 2 is the 
usual special relativity factor. In the notation of Sec. [TT1 
v cos 6 = v lOB . 

The Doppler shift effect can also be understood in two 
additional, complementary ways. The LISA response 
function naturally contains a Doppler shift term in the 
phase, due to LISA'S motion about the Solar System 
barycentcr. The Doppler shift discussed above is iden- 
tical to this, but now it is the GW source that moves, 
as opposed to the detector. Similarly, one could con- 
sider first the GW phase emitted in the EMRI's center 
of mass, and then map this to that observed in the So- 
lar System by shifting the phase's time-dependence by 
dt, corresponding to the light travel time along the line 
of sight between the center of mass of the EMRI system 
and the SMBH-SMBH system. From this perspective, 
the maximum phase shift that could accumulate in one 
year is simply the product of the EMRI orbital frequency 
and the light-crossing time of the projection of one year 
of SMBH binary evolution along the line of sight. 

We can then easily integrate Eq. (|A5j) , assuming a con- 
stant to, to recover the A$ GW computed in Scc|TTJ Setting 
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l = tt/2 = 5, we find 



w / v lo .(t)dt. 



N fM Se 



2 \ M T 



(A6) 

Notice that by choosing 6 = tt/2, there is no leading- 
order, unobservable constant velocity drift term. In the 



second line, we have Taylor expanded about w Ncwt t = 
and used that $ G w,Tot = N = ujT, where T is the time 
of integration, and that $ G w.Tot = N, where N is the to- 
tal number of radians in the non-accelerating waveform. 
Notice that this is the same A$ GW correction described 
in Sec. ED 
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